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Abstract 



Computational studies of basic models of strongly-correlated electron sys- 
tems can provide guidance in the search for new materials as well as insight 
into the physical mechanisms responsible for their properties. Here, we illus- 
trate this by examining what numerical computations have shown us about 
the 2-leg Hubbard ladder. 
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This workshop in honor of Zachary Fisk's 60 th birthday is to explore the future of ma- 
terials physics, an area to which he has contributed so much. It addresses the challenges 
and opportunities posed by complex materials. Here, in this session on computation, I want 
to focus on the problems of strongly interacting electronic systems, specifically the problem 
of the 2-leg Hubbard ladder. While this is a special case, I believe it illustrates the way in 
which numerical studies of simple models can help in the search for new materials, such as 
the cuprate ladder materials @,@], as well as how such calculations can provide insight into 
more general problems such as that posed by the high T c -layered cuprates. It also represents 
an example of a case in which the interest in these materials arose from numerical calcula- 
tions, which suggested that an undoped ladder would have a spin gap and a doped system 
would have power law pairing correlations ||. 

The Hamiltonian for the 2-leg Hubbard ladder is 

H = ~t ( C !+l,As C i,As + c[ Xs C i+1M ) ~t ± J2 (4,U<k,2s + c{ 2s C iM ) +Uj^ n i,X] n iM i 1 ) 
i,X,s i,s i,X 

Here c\ Xs creates an electron of spin s at rung % on leg A = 1 or 2. The hopping along a leg 
is t and across a rung is t±. There is an onsite Coulomb interaction U and the filling is set 
by a chemical potential in the Monte Carlo simulations and by a fixed electron number in 
the density matrix renormalization group (DMRG) calculations. 

The general properties of the 2-leg ladder M , which have been found from both numerical 
|5||| and analytic studies can be summarized as follows. At half-filling, the ground 

state has both a charge gap A c and a spin gap A s . It is an insulator with short-range 
antiferromagnetic correlations and spin gapped S = 1 magnon excitations. As t±/t increases 
the ground state of the half-filled ladder adiabatically crosses over from a spin-gapped Mott 
insulator to a band insulator || and A s /A c goes to 1. For small values of U, this crossover 
to the band insulator occurs when t±/t approaches 2. This corresponds to the value at 
which the bottom of the antibonding band rises above the top of the bonding band and the 
non-interacting, half-filled system becomes a Bloch insulator. When the spin gapped Mott 
insulator is doped so that the site occupation n = 1 — x, the system goes into a Luther- 
Emery phase [nj with power law pair field and 4/c^-CDW correlations. The sign of the 



rung-rung pair field correlations are positive while the rung-leg correlations are negative, 
implying "<i-wave-like" pairing. In the Luther-Emery phase, the power law decay of the 
pair field £~ 1 I 2k is conjugate to the £~ 2k decay of the 4/c^-CDW correlations. As one goes 
towards the insulating state, Umklapp processes play an important role |TT| and Schulz || 



has shown that k — > 1 and the pairing correlations are dominant. Finally, we note that the 
superfluid density vanishes as the doping x goes to 0. 

These properties make the 2-leg Hubbard ladder an interesting laboratory for studying 
the effect of various parameters on the pairing. In Figure 1, a DMRG calculation || of the 
rung-rung pair field correlation function 

D(£) = (A i+£ At) (2) 



is shown for a 2 x 32 ladder with a hole doping of x = 0.125, U/t = 8, and different values 
of t±/t. Here 

= ( c i,lT C i,2| ~~ C !,1| C !,2t) (3) 



2 



creates a singlet pair on the i th rung. As t±/t increases, the pairing correlations initially 
increase. However, as shown for t±/t = 2, when t±/t becomes too large, we find that the 
pairing correlations are suppressed. A closer look at this behavior is shown in Fig. 2a. Here 
an average strength of the rung-rung pairing correlations 

D = \Y.D{t) (4) 

is shown versus t±/t for U/t = 8 at several dopings. In Fig. 2b, D versus tj_/t is shown 
for various values of U/t at a doping x = .0625. Separate calculations of the pair gap show 
that it has the same behavior as D as a function of t±/t and x. The strength of the pair 
binding increases as x decreases and the Mott insulating state is approached while, as noted 
previously, the superfluid density decreases as x goes to zero. The pairing also depends 
upon U, and as shown in Fig. 2b, peaks for U of order the bandwidth. For values U large 
compared with the bandwidth, the effective exchange interaction varies as U^ 1 and both the 
insulating spin gap and the pairing strength decrease as U increases. 

The peak in D as a function of t±/t can be understood in terms of the behavior of 
the single-particle spectral weight A(p,u). Results for A(p,u) obtained from a maximum 
entropy analytic continuation of Monte Carlo data for a 2- leg ladder at T = .25t with 

U/t = 4, t±/t = 1.5, and (n) = 0.875 are shown in Fig. 3. Here, the solid curves are for the 
bonding band (p y = 0) and the dotted curves are for the antibonding band (p y = n). One 
sees that for this value of t±/t, the bonding band has spectral weight near the fermi level 
for p ~ (it, 0) while the antibonding band has its spectral weight near the fermi level for 
p ~ (0, 7r). These Fermi points can be connected by scatterings involving a large momentum 
transfer where, as we will see, the effective pairing interaction V is strong. Furthermore, the 
peak in both the bonding and antibonding spectral weight is seen to disperse very slowly 
near the bonding and antibonding fermi points leading to a large density of states at the 
fermi energy. For a value of t±/t ~ 1.6 which is slightly larger, the dispersion becomes even 
flatter and D peaks. At still larger values of t±/t, the antibonding band pulls away from 
the Fermi energy leaving only the bonding band with spectral weight near the Fermi level. 
In this case, the doped system behaves like a one-band Luttinger liquid and the pairing 
correlations vanish. 

These computations suggest that an array of weakly coupled ladders should exhibit 
superconductivity and that T c can be enhanced by adjusting the ratio of the effective rung- 
to-leg hopping as well as optimally doping the system. Additional calculations |H| have 



shown that the pairing can be further enhanced if an exchange interaction J±Sn ■ S 1 ^ is 
added. In a CuOi model, such an additional exchange term can be seen to arise from 
intermediate states in which the two exchanged electrons virtually occupy the oxygen site 
rather than one of the Cu sites. Chemically, one might try to enhance the rung hopping and 
the rung exchange by substitutions off the ladder plane that change the Madelung energy 
of the rung O. A more radical approach would be to replace the rung O by another element 
such as S. Clearly, the 2-leg ladder materials represent an important area of new materials 
that merit further study. 

Beyond this, however, computational studies of the 2-leg ladder have provided insight into 
the more general problem of the high T c cuprates. Using Monte Carlo simulations [|I~2] 



one 



can extract the effective pairing interaction V(q,u m ). While these calculations for a doped 
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ladder were restricted to temperature T > 0.25t because of the "fermion sign" problem, they 
showed that the momentum, Matsubara frequency and temperature dependence of V(q, uj m ) 
was remarkably similar to that of the spin susceptibility 

Xiq^m) = ^ [ dre-^J2e^(MUr)Mm) (5) 

iv JO q 

Here Mf = {n^ — n^) /2. A comparison of V(q,u m = 0) versus q with oj m = 0) for 
various temperatures is shown in Figures 4a and b. As the temperature decreases and the 
short-range antiferromagnetic correlations evolve in x(<7)> the effective interaction grows in 
strength and a broad peak appears at large momentum transfers. It is this interaction which 
leads to the strong scattering of pairs from regions near the bonding fermi points (pFb, 0) and 
(—PFb, 0) to regions near the antibonding fermi points (pFa, 7r) and (—pFa, —n) and hence 
to the pair binding. 

It is also useful to consider a strong-coupling local picture of the pairing. Figure 5 shows 
the results of a DMRG calculation of the amplitude for removing 2 electrons W& 



(N-2\(c^ Cjl -c iiCjl )\N) (6) 

from a half-filled 2-leg ladder. The phase of this amplitude has been chosen so that the rung 
amplitude is positive. Here one sees a clear "(4.2 _ y 2 -like" structure and from the size of the 
amplitude, one can judge that the pair is extended over about 10 sites. This amplitude is 
probing a property of the half-filled ground state and although the superfluid density for 
the undoped ladder vanishes, this spin-gapped Mott state contains latent d x 2_ y 2-\ike pairs. 
To illustrate this, consider a 2 x 2 segment of a ladder in which the 4-sites are numbered 
in a clockwise manner. In a ladder with short-range antiferromagnetic correlations, the 



half-filled, 4-site cluster has a ground state |L5 



^) = iY (A t 14 A5 3 -A t 12 Ay|0) (7) 



with A\j = cLct, — c]i c iT" ^ n * s su P er P os ition of valence bonds |L6| has d x 2_ y 2-\ike symme- 



try. The power law quasi-long-range-order superconducting pairing phase arises when these 
latent pairs are "freed" by the doping. 

In summary, the dependence of the pairing correlations on the ratio of the rung-to- 
leg hopping parameters t±/t, the strength of the Coulomb interaction U/t and the doping 
x provide insight into the search for new ladder materials that could exhibit enhanced 
superconductivity. The effective pairing interaction V(q,u m ) shown in Fig. 4a and the 
resonating valence bond structure shown in Fig. 5, illustrate the momentum space and real 
space nature of the pairing process for the 2-leg Hubbard ladder. The short coherence 
length and the relatively high energy scale set by the gap magnitude of the 2D layered high 
T c cuprates suggest that the basis pairing mechanism in these materials involves short-range 
physics. Thus we expect that what has been learned about the pairing process in the 2-leg 
ladder will also be relevant to the layered cuprates. 
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FIGURES 
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FIG. 1. The equal-time pair-field correlation function (AjAj) versus \i — j\ for (n) = 0.875, 
U/t = 8 and various values oit±/t. 
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FIG. 2. a) D versus i ± /t for J7/t = 8 at fillings (n) = 0.75, 0.875, and 0.9375. b) £> versus t±/t 
for various values of U/t at filling (n) = 0.9375. 




FIG. 3. Single-particle spectral weight A(p,uj) versus uo for t±/t = 1.5, T = 0.25t, U/t = 4, 
and (n) = 0.875. The solid curves denote the results for the bonding band (p y = 0) and the dotted 
curves denote the results for the antibonding band (p y = it). 
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FIG. 4. a) Momentum-dependence of the effective interaction V(q) for U = 4t, (n) = 0.875, 
and t± = 1.5i. Here V(q) is measured in units of t, q y = ir and V(q) is plotted as a function of 
q x . b) Momentum-dependence of the magnetic susceptibility x(q) for U = 4t, (n) = 0.875, and 
t± = 1.5t. Here, q y = it and x(q) i s plotted as a function of q x . 
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FIG. 5. Schematic drawing of the pair-wave function showing the values of the off-diagonal 
matrix element (^N — 2|(cj-fCjj — Cij_Cjf)| iV^ for removing a singlet pair between near-neighbor 
sites. 

-0.08 -0.08 

















+0.10 











-0.08 -0.08 



10 



